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MULTIMODE  SIMULATION  OF  HIGH  FREQUENCY  GYROTRONS 


I.  Introduction 

The  effort  to  push  gyrotron  oscillators  to  ever  higher  frequencies  and  power 
levels  has  led  to  the  employment  of  ever  higher  order  transverse  cavity  modes. 
Maximum  gyrotron  efficiency  is  generally  achieved  by  optimizing  the  magnetic 
detuning  for  a  single  mode.  This  becomes  progressively  more  difficult  to  do  as  the 
fractional  mode  spacing  shrinks,  since  1)  the  intended  mode  may  never  start 
oscillation,  or  2)  it  may  be  unstable  at  its  optimum  detuning  in  the  presence  of 
competition  from  nearby  modes,  or  3)  the  final  state  may  be  a  multimode  state.  In  this 
paper,  we  use  a  time-dependent  multimode  simulation  code  (Fliflet  et  al.  1991)  to 
examine  the  optimization  of  the  gyrotron  parameters  for  the  280  GHz  TE42,7-mode 
gyrotron  experiment  being  set  up  by  Kreischer  and  coworkers  at  the  Massachusetts 
Institute  of  Technology  (MIT),  in  the  presence  of  competition  from  nearby  TErn,l 
other  TE;n^  modes. 

A  classical  multimode  theory  of  lasers  was  obtained  by  Lamb  (1964)  and 
Sargent  et  al.  (1974).  The  full  nonlinear  system  of  coupled  equations  for  the 
interacting  modes  was  simplified  by  means  of  a  perturbation  expansion  in  powers  of 
the  electric  field  amplitudes.  This  theory  was  adapted  to  the  multimode  operation  of 
gyrotrons  by  Nusinovich  and  co-workers  (Moiseev  and  Nusinovich  1974,  Nusinovich 
1981)  who  also  made  use  of  a  perturbation  expansion  in  the  mode  amplitudes.  A  fully 
nonlinear,  time-dependent  multimode  theory  of  quasi-optical  gyrotrons  has  been 
developed  by  Bondeson  et  at  (1983).  The  stability  of  single-mode  gyrotron  operation 
with  respect  to  parasitic  oscillations  has  been  examined  by  Antonsen  et  al.  (1990), 
Levush  and  Antonsen  (1991),  Dumbrajs  et  al.  (1988),  and  Borie  and  Jddicke  (1988). 
Multimode  effects  in  gyrotron  have  also  been  studied  using  a  particle-in-cell  (PIC) 
code  by  Lin  and  co-workers  (Lin  et  al.  1988).  The  theoretical  basis  for  the  present 
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calculations  is  developed  in  Fliflet  et  al.  (1991).  A  fully  nonlinear  formulation  is 
desirable  v/ben  the  operating  current  is  many  times  higher  than  the  threshold  current. 
As  discussed  below,  there  is  a  difficulty  in  nonlinear  multimode  calculations  for 
configurations  involving  several  modes  with  different  azimuthal  and  radial  mode 
indices.  This  difficulty  is  associated  with  the  large  number  of  initial  phases  needed  to 
average  over  the  possible  relative  initial  phases  of  the  interacting  modes  which  in 
some  cases  can  increase  the  multimode  computation  time  by  a  factor  of  100-1000  over 
a  corresponding  single  mode  calculation.  Such  calculations  have  recently  become 
feasible  given  adequate  access  to  a  fast  computer  such  as  the  Cray. 

In  this  paper,  a  time-dependent  multimode  simulation  code  was  used  as  a  tool 
to  study  the  overall  mode  competition  problem,  including  competition  between  co-  and 
counterrotating  modes  for  the  same  m  and  I,  between  members  of  a  series  of 

modes  with  the  same  value  of  /,  and  between  modes  with  different  values 
of  m  and  1.  Here,  the  azimuthal  mode  index,  m,  is  generally  assumed  to  take  both 
positive  and  negative  integral  values.  However,  for  the  sake  of  clarity,  when  either 
+m  or  -m  is  explicitly  used,  it  is  assumed  that  m  is  a  positive  integer.  The  use  of  the 
plus  (minus)  sign  in  the  azimuthal  mode  index,  for  a  particular  TE±j„  g  mode,  refers  to 
the  sense  of  rotation  of  the  mode,  corresponding  to  co-  (counter-  )  rotation  of  the 
mode  with  respect  to  the  sense  of  electron  gyration  in  the  applied  axial  magnetic  field. 


11.  Outline  of  Theory 

The  total  transverse  electric  field  in  the  gyrotron  resonator  is  expanded  in  a 
series  of  transverse  modes  as  follows: 


N 

F=ya 
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where  an(t)  is  the  time  dependent  mode  amplitude,  h(z)  is  the  axial  profile  function,  e„ 
is  the  mode  transverse  vector  function,  coq  is  the  reference  wave  frequency,  and  V'n  is 
a  slow-time-scale  mode  phase  parameter.  The  free-running  oscillator  equations  for 
the  mode  amplitude  and  phase,  are  given  by: 


da, 
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where 


^n(0  =  d{oiQt)\^d\  h{z)e\  •  I 


(4) 


is  the  complex,  slow-time-scale  component  of  the  electron  beam  polarization  for  the 
rL  2 

mode  «,  W  =  dz\h{z)\  ,  cOnO  is  the  cold  cavity  eigenffequency  (taken  to  be  positive 
for  both  CO-  and  counterrotating  modes),  and  eo  is  the  free  space  permittivity.  To 
calculate  the  ac  current  density,  the  interaction  with  the  electron  beam  is  treated  in  the 
single-particle  approximation.  The  general  time-dependent  problem  can  be  simplified 
by  using  the  fact  that  the  characteristic  risetime  of  fields  in  the  resonator  is  much 
longer  than  the  electron  transit  time  in  the  cavity  as  well  as  the  wave  period.  In  this 
case,  one  can  use  a  quasi-steady-state  approximation,  in  which  the  electron 
trajectories  are  calculated  for  rf  modes  with  fixed  amplitude  and  linearized  phase. 

The  slow-time-scale  nonlinear  equations  of  motion  for  an  electron  in  a  thin 
annular  beam  immersed  in  a  tapered  magnetic  field  and  interacting  at  a  particular 
harmonic  with  several  circularly-polarized  TE  modes  have  been  given  previously  by 
Fliflci  (1986)  and  Fliflet  et  al.  (1982),  and  are  given  by: 


dz 


dz 
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where  u^  =  yVtl  c  is  the  normalized  transverse  momentum  amplitude,  u^  =  yo2l  c  is 
the  normalized  axial  momentum  and  Vz  denote  the  electron  transverse  and  axial 
velocities,  respectively), 

A={aQ-sniyQ)zlfi2  +  (OQtQi-s(p-{mz-s)eQ  (8) 

gives  the  slow  variation  in  the  transverse  momentum  azimuthal  phase  relative  to  the 
reference  wave  phase, 

¥n  =  ¥nM  +  ^A  (9) 

ToPzO 

is  a  phase  correction  corresponding  to  the  difference  between  the  mode  frequency  and 
the  reference  frequency,  rrir  is  the  azimuthal  index  of  a  particular  mode  called  the 
reference  mode,  c,  and  t=  (OQt.  Other  parameters  include  s,  the  harmonic 

number,  y  (yo),  the  (initial)  relativistic  mass  ratio  which  is  given  by: 

y  =  (l  +  u?+n2jl/2^  (10) 

kt ,  the  mode  transverse  wave  number,  ,  the  Larmor  radius  of  the  orbit,  7j(  Js),  the 
(derivative  of)  a  regular  Bessel  function,  Qz .  die  nonrelativistic  cyclotron  frequency 
for  the  z-component  of  the  magnetic  field,  Q,  the  value  of  i2z  at  the  resonator 
entrance,  toi,  the  time  the  electron  enters  the  cavity,  and  fn  ,  the  normalized  mode 

itn  A 

amplitude  with  azimuthal  dependence  e  "  .  The  mode  amplitude  for  >  0  is 
normalized  according  to: 
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(11) 
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The  sign  of  determines  the  direction  of  mode  rotation,  as  discussed  previously. 

The  normalized  wave  amplitude  for  <  0  is  given  by  Eq.  (11),  with  m„  replaced  by 
Im^l  and  replaced  by  (-l)^-f|^  Quantities  with  an  overbar 

have  been  normalized  according  to:  z  =  z  f  ry,,,  7^^  =  ri^  I  r^,  Q  =  Q  r^y/c ,  aQ  =  (OQ ry^/c , 
and  =  k^t  ry, .  Rq  and  0[i  denote  the  electron  orbit  guiding  center  radius  and 
azimuthal  angle,  lei  is  the  electron  charge,  wq  is  the  electron  rest  mass,  m„  is  the 
azimuthal  index  for  the  mode  n,  1*  a  zero  of  Jm\  f-n  is  the  radial  mode  index, 

Th,  is  an  arbitrary  normalization  radius,  and  (p  gives  the  slow  variation  in  the 
transverse  momentum  phase  relative  to  the  cyclotron  motion.  The  transverse 
TE-mode  normalization  coefficient  is  given  by: 


The  ac  current  density  is  obtained  by  integrating  Eqs.(5)-(7)  for  an  appropriate  set  of 
initial  conditions  at  the  cavity  input,  z=0.  For  a  thin  annular  beam,  the  transverse  ac 
current  density  is  given  by: 

(13) 

Using  the  prescription  developed  in  previously  in  Fliflet  et  al.  (1991),  the  mode 
amplitude  and  phase  can  be  rewritten  as: 


dr 
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where  (  Qq  denotes  the  average  with  respect  to  the  initial  momentum  phases  and 
guiding  center  angles  of  the  electrons,  and  L  is  the  normalized  interaction  length.  The 
normalized  current  is  given  by: 


7  _ 

- - 7 - .■■■■- — - - - —Iq  (15) 

moc^(Oo  n(l  -  mi  /  )W 

where  Iq  is  the  dc  beam  current,  and  the  frce-space  impedance  Zq=311  ohms.  The 
numerical  calculations  can  be  simplified  by  noting  that  the  argument  of  the  Bessel 
function  which  occurs  in  the  equations  of  motion  and  the  mode  amplitude  and 

phase  source  terms  can  be  expressed  as  -sv^l c,  and  therefore  these  Bessel 
functions  can  be  replaced  by  their  small  argument  expansion  with  little  loss  of 
accuracy. 

The  time-dependent  simulation  is  iititiated  by  assigning  a  small  initial 
amplitude  and  arbitrary  phase  to  a  set  of  modes  which  may  participate  in  the 
interaction.  The  corresponding  mduced  ac  current  density  is  obtained  by  integrating 
the  equations  of  motion  [Eqs.  (SM^)]  and  is  used  to  construct  the  source  terms  in 
Eqs.  (15)  and  (16).  Eqs.  (15)  and  (16)  can  then  be  integrated  for  a  single  time  step 
and  the  process  repeated.  The  initial  conditions  for  the  equations  of  motion  for  a  cold, 
phase-mixed  electron  beam  are:  w, (0)  =  w,o .  Mz(0)  =  ,  a  fixed  guiding-center 

radius  Rq\  and  yt(0)H/lo  and  ^  are  unifoiiuly  distributed  in  the  interval  [0,2;r].  The 
density  of  azimuthal  angles  must  be  high  enough  that  an  adequate  distribution  over 
this  interval  is  achieved  for  the  quantity  {m„  —mj)&Q  for  aU  values  of  mn-  In  the 
present  case  in  which  Imfll  and  Imnol  are  ~  40  and  can  have  the  same  or  opposite  sign, 
it  is  necessary  to  include  several  hundred  angles  to  preform  this  average  accurately  for 
all  the  interacting  modes.  The  interaction  efficiency  is  given  by: 

V  = ^ ^ (17) 

ro-i 
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and  the  output  power  in  the  mode  n  is  given  by; 


I 

JoriQC  , 
2ZQ\e\^ 


1- 


Qn-^rrin  ±s  ^ ) 


(18) 


for  modes  with  m„  =  ±l/n;,|,  respectively. 

For  comparison  with  other  work  it  is  useful  to  introduce  the  following 
normalized  parameters  which  are  often  used  in  gyrotron  analysis  (Danly  and  Temkin, 


1986): 


Fn  = 


2^~^s\ 


fn 


PzO  ^ 


fi  =  n 
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(20) 


(21) 


where  is  the  normalized  mode  amplitude,  /i  is  the  normalized  interaction  length,  and 
A  is  the  mode  detuning  parameter. 

III.  Conflguration  Modeled 

The  gyrotron  being  modeled  was  inspired  by  a  planned  MIT  280  GHz  gyrotron 
experiment.  We  were  provided  with  a  set  of  preliminary  calculations  for  that  device 
(Grimm  1991,  Kreischer  1991),  and  used  those  values  as  the  starting  point  for  our 
work.  The  gyrotron  is  designed  to  operate  in  the  TE42,7-mode  of  a  cylindrical  cavity 
with  wall  radius  /?h/=1.25  cm,  with  the  beam  placed  at  a  radius  rf,=0.747  cm,  in  order  to 
optimize  coupling  to  the  mode  co-rotating  with  the  electron  beam  (i.e.,  the  TE442,7 
mode,  in  our  notation).  Based  on  the  work  at  MIT,  the  optimum  magnetic  field  for  this 
gyrotron  was  5o=ll.l  T.  From  a  simulation  for  the  realistic  axial  profile  of  the  if 
electric  field,  MIT  derived  a  value  L/A=11.9,  where  L  is  the  normalized  interaction 
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length  (full  width  at  half  maximum)  in  the  cavity  for  an  equivalent  gaussian  profile,  and 
A  is  the  free-space  wavelength  of  the  radiation.  The  loaded  cavity  quality  factor  was 
Q=1 180.  The  beam  current  was  50  A  at  a  diode  voltage  of  90  kV.  However,  correcting 
for  6  kV  of  space  charge  depression  of  the  beam  voltage,  the  beam  energy  was  84 
keV.  The  anticipated  value  of  the  beam  a  is  1.6,  where  a  = 

In  general,  the  MIT  values  were  the  starting  point  of  our  study.  However,  the 
following  alterations  were  made.  First,  while  some  simulation  runs  have  been  canied 
out  at  the  indicated  value  of  L/A=l  1.9,  the  simulation  runs  reported  in  this  paper  were 
carried  out  at  the  slightly  shorter  value  L/A=8.26.  Second,  slightly  larger  magnetic 
fields  (smaller  detunings)  were  found  necessary  to  produce  stable  oscillation  either  in 
the  design  mode,  or  in  competing  modes  at  approximately  the  same  frequency,  as 
discussed  later  in  this  paper.  The  most  common  value  employed  in  the  simulations 
was  Bo=l  1.265  T. 

IV.  Results  of  the  simulations 

There  were  several  separate  issues  to  be  resolved  by  these  simulation 
studies,  and  a  large  parameter  space  to  be  investigated.  In  order  to  conserve 
simulation  time,  many  of  these  issues  were  tackled  separately,  in  limited  simulation 
runs  including  only  a  small  number  of  modes,  a  limited  voltage  ramp,  or  a  limited 
number  of  simulation  particles.  During  these  shorter  runs,  the  regime  of  stable 
operation  in  a  particular  final  state  was  established,  and  the  nature  of  the  mode 
competition  clarified.  Finally,  this  led  to  a  set  of  “long”  simulations,  which  started  at 
low  voltage  and  attempted  to  model  the  entire  sequence  of  modes,  from  the  first  that 
begins  oscillating  up  to  the  final  state. 

Because  of  the  cost  of  such  long  simulations,  compromises  were  necessary 
with  respect  to  the  various  simulation  constants.  An  important  limiution  of  this  smdy 
was  the  various  compromises  required  to  keep  the  running  time  of  the  simulation 
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within  available  resources.  There  are  a  number  of  simulation  constants  whose  scale 
size  will  affect  the  accuracy  of  the  simulation.  These  include  1 )  the  spabal  mesh  used 
to  evaluate  the  interaction  in  the  cavity,  2)  the  time  step  used  to  advance  the  fields,  3) 
the  slope  of  the  voltage  ramp  used  to  model  start-up,  4)  the  number  of  particle  phases 
included  at  each  azimuthal  location,  and  5)  the  number  of  azimuthal  beam  positions 
(particles)  included.  During  the  shorter  runs,  an  attempt  was  made  to  establish 
“experimentally”  the  minimum  values  of  various  simulation  constants  required  to 
pioduce  a  simulation  that  was  “qualitatively”  in  agreement  with  the  result  of 
simulations  employing  higher  numbers  of  particles,  more  mesh  poi’^^s,  etc.  Based  on  a 
number  of  computational  tests,  these  compromises  should  not  affect  the  resulting  final 
principal  mode  in  the  long  simulations.  However,  the  mode  purity  of  the  final  state 
was  separately  addressed  by  using  finer  simulation  constants  in  more  limited 
simulation  runs. 

Additional  limitations  of  this  study  include  the  use  of  a  beam  with  no  guiding 
center  or  a  spread,  and  the  (arbitrary)  modeling  of  the  beam  current  and  a  during  the 
rise  of  the  voltage  pulse.  The  presence  of  guiding  center  spread  would  further 
complicate  the  mode  competition,  since  additional  series  of  modes  would  strongly 
couple  to  the  beam,  and  the  relative  coupling  of  the  various  modes  would  change. 

Such  spread  could  be  included  in  the  model  in  a  straightforward  way,  by  averaging  the 
coupling  coefficient  for  each  mode  over  the  radial  spread — however,  this  was  not  done, 
'fhe  presence  of  a  spread  would  slightly  modify  the  growth  rates  of  the  various  modes, 
and  also  have  a  small  effect  on  the  saturated  efficiency.  Including  a  spread  would 
increase  the  particles  required  in  each  simulation  by  an  order  of  magnitude  (In  the 
experiment,  it  might  also  limit  the  maximum  average  a  achievable.)  This  was  not 
done.  The  current  and  a  modeling  are  discussed  in  the  next  section. 
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N.l  Results  of  the  preliminary  simulations 

Before  examining  the  simulation  results,  it  is  useful  to  examine  calculations  of 
the  threshold  for  oscillation  (starting  current)  for  various  modes.  Starting  current 
calculations  were  carried  out  using  an  analytic  model  that  employs  a  sinusoidal 
approximation  to  the  axial  profile  function  of  the  mode.  In  order  to  match  this  to  the 
gaussian  model  used  by  the  time-dependent  code,  the  length  and  cavity  Q-factor  were 
adjusted  from  the  nominal  values  used  in  the  time-dependent  simulation  to  cause  the 
minimum  starting  current,  and  its  corresponding  voltage,  to  agree  with  that  for  a 
gaussian  profile  for  the  TE+42,7  mode.  The  starting  current  calculations  were  carried 
out  for  all  first  harmonic  modes  with  |  m  I  ^0  and  where  the  third 

subscript  “1”  refers  to  the  lowest  order  axial  mode. 

Starting  current  calculations  as  a  function  of  voltage  (beam  kinetic  energy)  for 
Bo=l  1-265  T  and  a  normalized  radius  rhlR^=Q.6,  corresponding  to  a  beam  radius 
r^^OJ^l  cm,  arc  shown  in  Fig.  1.  In  these  calculations,  the  beam  a  is  proportional  to 
beam  voltage,  with  05=1.6  at  84  keV,  while  the  beam  current  is  assumed  to  be  50  A, 
independent  of  voltage.  Also,  the  Q  of  each  mode  is  assumed  to  scale  as  the  square  of 
the  mode  frequency,  with  Q=1 180  for  the  TE42J  mode.  This  modeling  of  current  start¬ 
up  is  somewhat  realistic,  since  the  electron  beam  diode  runs  temperature  rather  than 
space-charge  limited,  and  the  current  is  known  experimentally  to  turn  on  much  more 
quickly  than  the  voltage.  The  assumption  that  a  is  proportional  to  the  beam  kinetic 
energy  is  reasonable,  given  that  the  initial  a  in  the  diode  is  produced  by  the  initial 
accelerati(;n  of  the  electrons  by  the  electric  field  with  a  component  normal  to  the 
magnetic  field.  However,  no  attempt  has  been  made  to  model  the  current  and  beam  a 
more  realistically  through  such  means  as  particle  simulations.  These  same  scaling 
laws  were  used  in  the  time-dependent  simulations. 

For  this  TE+42,7-mode  design,  the  design  mode  has  the  lowest  starting  current 
(and  highest  growth  rate)  for  a  limited  range  of  voltages  in  the  vicinity  of  70  keV.  At 
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lower  voltages,  other  TE+mj  and  TE_ff,,8  modes  have  the  lowest  starting  currents  (and 
highest  growth  rates).  Based  on  a  slow  voltage  risetime,  the  system  should  be 
expected  to  tune  through  a  sequence  of  higher  frequency  modes  before  reaching 
conditions  under  which  the  TE+42,7  mode  can  start  oscillation.  When  the  beam  line 
enters  the  start  oscillation  curve  for  this  mode,  oscillation  can  build  up  provided  that  it 
is  not  suppressed  by  the  rf  electric  fields  of  another  mode.  In  a  simple  case,  in  which 
the  only  relevant  modes  were  a  sequence  of  TE+^j  modes,  the  TE442,7  mode  would 
be  expected  to  grow  to  saturation  when  the  system  demnes  out  of  the  TE+43  7  mode. 

It  is  noteworthy  that  the  50  A  beam  current  is  greater  than  10  times  the  starting 
current  at  70  keV  for  the  TE442.7  mode,  and  in  fact  exceeds  the  starting  current  for 
more  than  15  other  modes.  However,  the  linear  growth  rate  of  each  mode  is  inversely 
proportional  to  the  starting  current,  so  that  for  a  system  starting  from  noise  at  70  keV, 
the  TE+42,7  mode  should  oscillate. 

Accordingly,  the  first  simulations  employed  just  the  TE442,7  mode,  as  in  the 
MIT  design,  and  determined  its  optimum  detuning  (i.e.,  magnetic  field),  and  its  final 
power  and  efficiency.  The  next  simulations  employed  just  the  TE+^jj  series  of  modes 
and  determined  the  effect  of  mode  competition  and  multimode  effects.  Mode 
competition  limited  the  maximum  detuning  for  the  TE442,7  mode  to  the  point  at  which 
the  final  state  became  the  TE+41J,  causing  some  reduction  in  efficiency,  while 
multimode  effects  limited  the  purity  of  the  final  state,  as  significant  TE+41  j  and 
TE+43  7-mode  satellites  appeared  in  the  final  TE+42,7  state. 

Each  of  the  TE±,„,/  modes  has  a  mode  accompanying  it  at  the  same 

frequency,  but  with  the  opposite  sense  of  rotation  and,  in  general,  a  different  coupling 
to  the  beam.  The  next  issue  to  be  examined  was  the  effect  of  including  the  TE_;„,7 
series  of  modes  in  the  simulations.  (This  required  a  large  increase  in  the  number  of 
simulation  particles,  as  discussing  in  Section  2.)  These  oppositely  rotating  modes 
compete  with  the  TE+n,l  modes,  but  have  a  lower  growth  rate  for  the  0.747 -cm  beam 
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radius  used  in  the  simulations.  It  was  found  that  the  final  state  always  consisted  of 
modes  with  just  one  sense  of  rotation,  as  the  stronger  mode  quickly  suppressed  the 
weaker  one  of  opposite  rotation.  Thus,  the  mode  competition  between  a  TE+m,/  mode 
and  a  TE+^i  /  mode  does  not  lead  to  a  multimode  state  (i.e.,  a  standing  mode).  (In 
fact,  prestarting  the  simulation  with  a  large  signal  in  the  sense  of  rotation  with  the 
lower  growth  rate  can  lead  to  apparently  stable  oscillation  whose  sense  of  rotation  is 
opposite  to  that  which  occurs  due  to  ordinary  start-up  from  noise.)  However,  during 
the  mode  hopping  that  occurs  during  the  voltage  ramp,  short  intervals  of  growth  of  the 
TE_;„  7  modes  could  be  observed,  before  the  corresponding  (and  faster  growing) 

TE+;„  7  modes  grew  large  enough  to  suppress  them. 

The  only  effect  of  keeping  the  TE_;„  j  modes  in  the  simulations  was  an 
apparent  slight  reduction  in  the  initial  growth  rate  of  the  TE+;„j  modes.  In  a  situation 
involving  strong  competition  between  modes  with  almost  identical  growth  rates  and 
different  values  of  m  and  i,  however,  the  presence  of  these  TE_;„  j  modes  can 
potentially  affect  the  final  state.  Simulations  including  only  the  t-7  series  of  modes, 
but  with  m=±40,  ±41,  ±42,  ±43,  ±44,  have  found  a  case  with  a  magnetic  field 
fio=l  1.2375  T  and  a  voltage  ramp  of  1.8  kV/ns,  in  which  the  system  will  skip  from  the 
TE443  j  to  the  TE441J  due  to  the  inclusion  of  oppositely  rotating  modes,  while  it  ends 
up  in  the  TE+42,7  if  only  the  plus  modes  are  used  in  the  simulation.  This  appears  to  be 
a  risetime  effect  in  the  present  case.  However,  it  is  easy  to  imagine  cases  with  more 
closely  spaced  modes  in  which  the  TE^l  /  mode  will  be  skipp>ed  over  in  favor  of  the 
TEm-2/  mode,  even  over  very  slow  voltage  ramps. 

The  next  issue  to  be  resolved  was  the  affect  of  including  TE;„^  modes  with 
based  on  Fig.  1,  the  most  troubling  modes  were  in  the  TE_^  8  series.  These 
modes  in  general  occur  at  slightly  higher  frequencies  than  the  nearest  TE+^  7  modes, 
but  have  comparable  coupling  to  the  electron  beam  at  the  design  radius,  and  thus  have 
a  temporal  advantage  during  the  rise  of  the  voltage  pulse.  In  fact,  for  simulation  runs 
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beginning  at  much  lower  voltages,  this  series  of  modes  invariably  suppressed  the 
TE+ot  j  modes,  even  though  the  TE+42,7  mode  has  the  lowest  start  oscillation  current 
in  the  vicinity  of  the  final  voltage.  The  reason  for  this  will  be  discussed  in  this  next 
section. 

rv.2  Results  of  the  long  simulations 

The  culmination  of  a  sequence  of  short  runs  to  investigate  parameter  space  and 
to  understand  the  effect  of  varying  the  various  simulation  constants  was  the  first  long 
computer  run  (~9  hours  of  running  time  on  a  Cray  X-MP)  that  attempted  to  model  the 
entire  time  sequence  during  the  voltage  ramp,  from  start-up  of  the  first  mode,  to  the 
final  state  at  full  voltage  (84  keV).  This  simulation  is  shown  in  Fig.  2.  The  beam 
radius  was  chosen  to  correspond  to  the  MIT  value  irblR^=0.6)  that  optimizes  coupling 
to  the  TE42,7  mode,  even  though  the  result  is  operation  in  the  TE_39,8  mode. 

However,  the  magnetic  field  was  chosen  as  Bo=l  1-265  T  to  produce  stable  oscillation 
in  the  TE_39,8  mode  at  84  keV.  Based  on  an  examination  of  Fig.  1,  and  the  results  of 
more  limited  simulation  runs,  only  the  TE+m,l  and  TE_ff,,8  modes  were  included  in  this 
simulation,  with  the  oppositely  rotating  modes  in  the  1=1  and  /=8  series  excluded.  A 
tiO-ns  voltage  ramp  was  employed,  beginning  at  30  keV,  followed  by  a  30-ns  voltage 
flat-top  at  84  keV.  The  simulation  constants  were  chosen  to  minimize  the  length  of 
the  simulation  without  affecting  the  final  results  in  a  significant  way.  (However,  the 
mode  purity  of  the  final  state,  and  the  final  efficiency,  were  the  subject  of  more  limited 
runs  with  much  finer  simulation  constants  that  confirmed  the  efficiency  and  multimode 
content  of  the  final  state.) 

Figure  2  shows  that  the  first  mode  to  reach  high  power  is  the  TE_43,8  mode. 
(While  Fig.  1  suggests  that  the  TE_44,8  mode  should  oscillate  at  ~26-28  keV,  this 
was  not  seen  in  the  preliminary  simulations,  and  more  accurate  starting  current 
calculations  using  the  same  gaussian  profile  used  in  the  simulations,  suggest  that  the 
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actual  starting  current  for  this  mode  is  ^  50  A.)  Shortly  after  the  start  of  the  TE^3  g 
mode,  the  TE+46  j  mode  begins  to  grow,  but  is  quickly  suppressed  by  the  large 
amplitude  TE_43,8-mode  oscillation.  However,  two  satellite  modes,  the  TE^2,8  mode 
and  the  TE  44^8  mode,  grow  up  at  about  10  ns.  Subsequently,  the  TE_43,8  mode  is 
detuned  by  the  voltage  ramp,  and  begins  to  fall  off.  The  higher  frt  rency  TE_44,8 
satellite  mode  also  falls  off.  However,  the  lower  frequency  satellite,  the  TE_42,8  mode, 
begins  to  grow,  saturates  at  ~16  ns,  and  renuiins  at  high  power  till  ~24  ns.  During 
this  period,  the  previous  oscillating  mode,  the  TE^3  g  mode,  has  become  a  satellite 
mode,  and  a  new  satellite,  the  TE_4i,8  mode,  has  appeared.  This  mode  progression 
continues  through  the  voltage  ramp  until  the  final  TE_393  mode  oscillates.  Its 
satellites  are  the  TE_4o,8  mode  and  the  TE_38_8  mode.  Figure  3  shows  the  overall 
multimode  efficiency  during  the  mode  sequence. 

Figure  2  shows  mode  evolution  within  the  family  of  TE_ffi,8  modes,  with  no 
serious  competition  from  the  TE+^j  modes.  Each  of  the  intermediate  states,  from  the 
first  mode  to  stan  oscillation,  the  TE_43,8  mode,  to  the  final  TE_393  mode,  is  a 
multimode  state,  with  its  m'=m±\  mode  satellites  present  at  levels  >1%.  However, 
close  examination  of  the  simulation  shows  that  active  mode  competition  takes  place 
only  during  each  transition  from  a  particular  TE_ot,8  mode  to  the  lower  frequency 
TE_^1,8  mode.  This  competition  involves  both  TE_^,8  modes  and  TE+^jj  modes. 

This  can  be  seen  in  Fig.  4,  which  shows  only  the  TE+;„  7  modes  from  the  simulation. 
With  the  exception  of  the  TE+46^7  mode,  the  remaining  TE+mj  modes  show  growth 
spurts  only  during  the  transition  between  oscillation  in  adjacent  TE_;„^8  modes.  They 
are  mmed  loose  to  grow  as  each  TE_;„,8  mode  detunes  out  of  its  high  power 
oscillation,  but  are  quickly  suppressed  as  the  next  TE_^+l^8  mode  grows  from  a 
satellite  mode  to  high  power  oscillation.  In  the  final  state,  the  power  in  the  main 
TE_39,g  mode  is  -1.2  MW,  while  the  TE_4o,8  satellite  is  at  45  kW  and  the  TE_38,8 
satellite  is  at  30  kW.  Thus,  the  overall  mode  purity  is  about  94%. 
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Figxire  5  shows  the  result  of  an  identical  simulation  employing  only  the  TE+m,l 
modes.  As  expected,  the  system  proceeds  in  a  sequence  of  TE+m,?  modes,  and  the 
final  state  is  the  TE+42,7  mode,  with  TE441  j  and  TE+43  j  satellite  modes. 

5.  Discussion  of  the  results 

It  was  found  that  for  systems  of  the  order  of  length  of  the  proposed  MIT  design 
(or  somewhat  shorter),  the  final  state  (i.e.,  during  the  voltage  flat-top)  was  a 
multimode  state,  with  the  m±l  satellites  of  the  oscillating  TE;„,/  mode  present  at  the 
~2-5%  power  level.  This  has  some  consequence  with  respect  to  the  utility  of  the  tinal 
state  for  particular  applications.  In  addition,  the  intermediate  states  during  the 
start-up  voltage  ramp  were  multimode  states,  an  effect  that  plays  a  strong  role  in 
determining  the  particular  mode  which  will  dominate  the  final  state.  The  multimode 
complexion  of  these  oscillating  states  seems  to  be  a  function  of  the  mode  density.  To 
test  this  hypothesis,  a  simulation  run  was  carried  out  employing  only  alternate  TE_ff,,g 
modes.  In  this  case,  the  satellite  modes  have  twice  the  frequency  separation.  The 
maximum  satellite  level  seen  was  <0.1%. 

Mode  competition  is  critical  during  the  voltage  risetime.  Generally  speaking, 
due  to  the  presence  of  a  voltage  ramp  at  the  leading  edge  of  the  voltage  pulse,  the 
operating  mode  is  not  the  first  mode  to  reach  high  power,  but  the  last  mode  to  reach 
high  power,  during  the  voltage  ramp.  More  specifically,  the  system  passes  through  a 
sequence  of  high  power  modes.  (In  the  simulations,  these  intermediate  modes  last 
only  ~10  ns.  However,  in  the  much  slower  experimental  voltage  ramp,  each  mode 
might  last  hundreds  of  ns.)  The  same  type  of  mode  sequence  was  invariably  observed 
during  the  voltage  ramp,  as  follows.  The  first  mode  to  start  oscillation  simply  started 
a  sequence  of  ~5  high  power  modes,  at  progressively  lower  frequency,  as  each  mode 
was  detuned  with  increasing  voltage,  leading  to  the  final  mode.  (Final  here  must  be 
qualified  by  noting  that  we  have  modeled  only  a  few  tens  of  nanoseconds  of  the 
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voltage  flat-top,  and  the  apparent  stability  during  this  time  interval  may  not  persist  for 
the  microseconds  of  the  experimental  voltage  pulse.)  We  have  found  that  each 
intermediate  state  ir  a  multimode  state,  since  when  a  TEfnj  mode  is  large,  it  will  have 
satellites  at  TE^ij.  This  is  a  result  of  a  "well-understood"  process  (Levush  and 
Antonsen  1990)  whereby  the  modes  grow  up  due  to  a  resonance  condition 

involving  However,  TE^'^t  modes  in  other  radial  series 

with  mVm±l,  do  not  become  satellites,  and  in  fact  are  found  in  the  simulations 
to  be  stronglv  suppressed.  Since  the  satellite  modes  are  from  the  same  set  of  radial 
modes  (e.g.,  TE+„j  or  modes),  the  lower  frequency  TE+^^+l^/  satellite  mode 

will  have  both  gain  and  higher  initial  amplitude  than  other  competing  modes,  just  as 
the  TE+^^  mode  detunes  out  of  its  nonlinear  operation  due  to  the  effect  of  the  positive 
voltage  ramp. 

In  fact,  the  simulations  indicate  the  situation  is  still  more  strongly  biased  in 
favor  of  the  lower  frequency  satellite  mode.  That  is,  the  TE+wfi,/  mode,  which  begins 
at  ^  1%  of  the  power  of  the  main  mode  (or  ^  10%  of  its  rf  electric  field),  will  begin  to 
grow  and  ultimately  to  suppress  the  TE+^i^  mode,  due  to  the  voltage  ramp,  so  that  the 
TE±fn^l  mode  acmally  turns  off  at  a  slightly  lower  voltage  than  it  would  in  a  single 
mode  run.  Thus,  no  other  modes  have  much  likelihood  of  taking  over,  and  there  is  a 
strong  tendency  for  the  system  to  evolve  in  a  sequence  of  modes  with  constant  radial 
index  (and  constant  sense  of  rotation),  even  if  it  would  appear  that  at  full  voltage,  this 
set  of  modes  does  not  have  the  lowest  threshold  current  (highest  growth  rate). 

There  is  a  window  in  time  during  which  the  dominant  mode  has  begun  to  fall  off, 
and  the  lower  frequency  satellite  mode  to  grow,  during  which  other  modes  can  also 
compete.  However,  it  is  a  very  limited  window.  That  is,  while  the  satellite  is  growing 
by  only  -lOOx  in  power  to  saturation  (and  only  ~30x  to  the  point  at  which  nonlinear 
suppression  of  other  modes  begins),  another  mode  would  have  to  come  out  of  the 
noise  level  (which  was  arbitrarily  set  at  1  kV/cm  at  the  beam  for  each  mode),  and 
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grow  large  enough  to  suppress  the  satellite  mode.  Since  this  competition  is  based  on 
relative  growth  rates,  rather  than  the  voltage  risetime,  it  is  not  necessary  to  simulate 
a  ramp  as  slow  as  the  experimental  ramp  (several  microsecond  risetime)  to 
adequately  model  this  aspect  of  mode  competition. 

Therefore,  to  design  for  operation  in  the  TE442,7  mode,  for  instance,  it  is  not 
sufficient  to  ensure  that  this  mode  has  the  lowest  starting  current  in  the  vicinity  of  the 
final  voltage.  In  practice,  at  the  beam  radius  chosen  to  optimize  the  coupling  to  the 
TE+42,7  mode,  the  TE-^.s  modes  always  stan  first,  and  thus  the  final  state  is  instead 
the  TE_39^8,  even  though  this  mode  has  a  higher  starting  current  in  the  vicinity  of  the 
design  voltage.  Instead,  one  must  ensure  that  a  higher  frequency  TE+;„  7  mode,  such 
as  the  TE+45  7  mode,  is  the  first  mode  to  start  oscillation  during  the  rise  of  the  voltage 
ramp. 

To  test  this  supposition,  a  search  was  made  for  a  beam  radius  at  which  a 
TE+nj  mode  would  oscillate  first.  Figure  6  shows  the  starting  current  as  a  function  of 
voltage  for  the  case  of  1-265  T  and  rblRy^^.blA,  corresponding  to  a  beam  radius 
of  0.780  cm.  In  this  case,  the  TE+47  j  mode  should  be  the  first  mode  to  oscillate, 
followed  by  the  TE446,7.  (However,  there  is  a  small  interval  near  32  keV  where  the 
TE_43  g  mode  has  the  lowest  starting  current.)  Figure  7  shows  the  results  of  a 
time-dependent  simulation  run  for  this  beam  radius.  (The  remaining  parameters  are 
as  in  Fig.  2.)  As  expected,  the  TE+ff,,7-mode  series  dominates  in  a  sequence 
beginning  with  the  TE47,7  mode  and  ending  with  the  TE442,7  mode.  The  TE_43,8  mode 
shows  some  growth  at  ~10  ns,  but  is  quickly  suppressed  by  the  TE446,7  mode.  This 
simulation  run  included  the  TE_ff,j  modes  as  well,  since  the  TE_42,7  mode  actually  has 
a  lower  starting  current  than  the  TE442,7  mode.  However,  the  progression  is  a  series 
of  TE+ff,  j  modes,  as  predicted  from  the  previous  long  simulation,  and,  despite  its 
higher  linear  growth  rate,  the  TE_42,7  mode  is  completely  suppressed  in  the  final 
state,  because  the  TE.t42,7  mode  was  already  present  at  a  significant  level  as  the 
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system  detuned  with  voltage  from  the  TE+43,7  mode.  The  saturated  state  again 
appears  to  be  a  multimode  state,  with  substantial  amounts  of  the  TE+43J  and  TE+41  j 
modes  present. 

6.  Summary 

We  have  carried  out  a  simulation  study  of  mode  competition  in  highly 
overmoded  gyrotron  oscillators.  The  parameters  of  our  simulation  were  chosen  in  light 
of  a  280  GHz,  1  MW  gyrotron  design  from  MIT.  The  study  addressed  the  difficulty  of 
causing  the  system  to  oscillate  in  the  design  mode  in  the  presence  of  competition  from 
other  modes,  and  demonstrated  the  difficulty  of  achieving  stable  single-mode 
operation.  While  the  simulations  were  computationally  intensive,  and  could  not 
pursue  all  issues  relevant  to  long  pulse  operation,  they  were  used  to  study  numerous 
separate  cases,  and  established  trends  with  respect  to  the  problems  of  mode 
competition  during  the  start  of  the  voltage  pulse.  They  also  addressed  the  mode 
purity  of  the  final  state.  Specifically,  it  was  found  that  because  of  multimode  coupling, 
the  system  had  a  very  strong  tendency  to  remain  in  a  particular  series  of  TE^^  modes, 
with  m  of  a  particular  sign,  and  I  constant,  during  the  voltage  ramp.  Thus,  controlling 
the  first  mode  to  oscillate  can  permit  control  of  the  final  mode  to  oscillate,  even  in  the 
presence  of  substantial  competition  from  modes  with  different  values  of  1.  The 
efficiency  of  the  final  state  depends  strongly  on  the  final  detuning  (as  well  as  on  the 
presence  of  parasitic  modes).  In  general,  the  optimum  detuning  derived  for  a  single 
mode  was  unachievable  due  to  mode  competition.  Finally,  for  systems  of  the  length  of 
the  MIT  design  (or  somewhat  shorter),  the  final  state  appears  to  be  a  multimode 
state,  with  the  TE^tti^  modes  present  at  approximately  the  3-5%  level.  This  will 
degrade  the  useful  efficiency  of  the  gyrotron  for  many  applications.  In  addition,  the 
efficiency  is  somewhat  degraded  because  the  magnetic  field  detuning  that  optimizes 
efficiency  in  a  single  mode  run,  will  not  generally  result  in  operating  in  the  same  mode 
in  the  presence  of  mode  competition. 
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Figure  3.  Efficiency  from  the  simulation. 
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